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A possibility ol spontaneous magnetization in high density symmetric quark matter is 
investigated using the NJL type effective model of QCD by means of an effective potential 

■ with respect to an auxiliary field. It is shown that the quark ferromagnetic condensate has 
non-vanishing value at high baryon density due to the tensor-type four-point interaction 
between quarks. 

§1. Introduction 

O '. 

^ . One of recent interests about the field governed by the quantum chromodynam- 

Ch I ics (QCD) may be to understand the phase structure on the plane depicted by the 

temperature and the baryon chemical potential. This research area is stimulated 
CN , by the recent high energy heavy ion collision experiments, in which it is reported 

that the quark-gluon phase is realized.'^ Further, it is reported in the astrophys- 

■ ical area that there exist neutron stars with very large magnetic field which are 
, called magnetars.® In these extreme conditions at high density and/or at high 

temperature in the hadronic and/or the quark-gluon matter, there may exist var- 



(N 



X 



' ious phases such as hadron phase, quark-gluon phase, chiral symmetric or broken 

, phase,-- two flavor color-superconducting phase, color-flavor locked phas^ and so 

forth. The existence of a ferromagnetic phase has been previously considered by 
some authors .'^''^''^''^''^''^ In the investigation of the phase structure, especially 
in the system at finite density, various effective models of QCD are used because it 
is difficult to compute physical quantities at finite density directly from QCD, for 
a ■ example, from the lattice QCD simulation method. 

The Nambu-Jona-Lasinio (NJL) modeP^ as an effective model of QCD gives 
various useful informations for the system at finite temperature and density.'^ ''^''^ 
Especially, as for the chiral symmetry breaking and restoration, the NJL model has 
taught various properties until now. Further, the existence of the color supercon- 
ducting phase has also been investigated by using the NJL model. In this paper, 
with the extension of the original NJL model, we investigate a possibility of the spon- 
taneous spin polarization or magnetization in quark matter at very high density. 



typeset using VTpT^.ds (Ver.0.9) 



2 



Y. Tsue, J. da Providencia, C. Providencia and M. Yamamura 



Here, let us consider the NJL-type Lagrangian density which has chiral sym- 
metry, namely, the Lagrangian density is invariant under the su{2)l x sn(2)/j-chiral 
transformation. Within the four-point interaction between quarks, it is possible for 
a general form of the Lagrangian density with chiral symmetry to include scalar, 
vector and tensor interactions which is written as 

^ = -^kin + £s + + , 

Cs = -Gs [(^V)' + i^pij^ri^f] , 

Cv = -Gv [{i'l^Til^f + (V^Tst'^^V')'] , 

/:t = -Gt [(V^T^T'rV')' + (MsT^t'V')'] • (1-1) 

Here, ifj represents the quark field and r represents the isospin operator. The orig- 
inal NJL model Lagrangian density is £njl = >Ckin + ^5- Here, 'ipili and ipi'y^T'ip 
correspond to the sigma meson (a) and the pion (vf) filed, respectively, in terms of 
the bosonized meson fields. As is well known, at low energy, the chiral symmetry 
dynamically breaks down and the chiral condensate which is the vacuum expectation 
value of ipip becomes non-zero. Similarly to the original NJL model, the field ip'j^^Tip 
and -ip^^j^fip in Cy corresponds to the p meson and a meson fields, respectively, if 
these composite fields are treated by means of a bosonization. However, it is known 
that the physical a meson field is constructed from the original a meson field by 
mixing with vr and p mesons due to the chiral symmetry breaking with the chiral 
condensate (V'V')) which results to the mass difference between p and a mesons.'^ 
In addition to the above four-point interactions, further, from the structure of the 
Dirac gamma matrices, the tensor interaction part, Ct, should be introduced. Here, 
ipj^j'^Tip may correspond to the 02 meson field. In the nonrelativistic limit, this ten- 
sor meson exchange interaction by 02 meson leads to the tensor force, the spin-orbit 
force, the second spin-orbit force and so on.13 This term may also be understood 
in the framework of a Fierz transformation of a standard NJL interaction. 

In this paper, we consider a system at high density such as the interior of the 
neutron star, in which there may exist quark matter at high density, and/or the quark 
star .11^ At high baryon density, the chiral symmetry is restored and the dynamical 
quark mass is zero, namely, the chiral condensate {ipip) is equal to zero. In this case, 
the non-vanishing vacuum expectation values of V'V', ''pilbTip, ipj^r^ and ip^^^^rtp 
do not appear, that is, the condensates are zero for these composite field operators. 
Especially, as for the spin polarization originated from the axial vector interaction, 
it may be possible to realize the spin polarization {■ip^^'^'^Til))^^^ However, at high 
baryon density, where the chiral symmetry is restored and the quark mass is zero, 
the spin polarization does not appear.El' Therefore, it is regarded that these terms 
only cause the excitation modes around the vacuum which correspond to meson 
excitation modes of o", vr, p and a mesons in the low energy vacuum with a finite 
chiral condensate. Thus, we could safely discard the £5 and Cy when we only 
investigate the vacuum at high density 

In RefsJT8|l and IT9|) . by two of the present authors (J. da P. and C. P.), a 
possibility of the ferromagnetic condensate in the quark matter at high density has 
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already been considered by means of the thermodynamical treatment of the NJL 
model with tensor interaction. In this paper, we reinvestigate the ferromagnetic 
condensate by a field theoretical approach with the aim of the extension to the finite 
temperature system. Namely, we derive the effective potential for the ferromagnetic 
condensate. 

This paper is organized as follows. In the next section, the effective potential for 
the quark ferromagnetic condensate is derived by using the auxiliary field method. 
In §3, the analytical form of the effective potential is given at zero temperature. 
In §4, the numerical results are given and it is shown that the quark ferromagnetic 
condensate is really realized for high density quark matter. Further, the implication 
to the magnetars is discussed. The last section is devoted to a summary. In Appendix 
A, the effective potential is reinvestigated in terms of another normalization. In 
Appendix B, the effect of the vacuum polarization is summarized briefly. 

§2. Effective potential approach to symmetric quark matter with 

tensor-type interaction 

In this section, we derive an effective potential with respect to the quark ferro- 
magnetic condensate. As is mentioned in §1, since we pay attention to the vacuum 
at high density quark matter, the chiral condensate and the excited modes may be 
discarded safely. Thus, let us start with the following Lagrangian density: 

jr = z^T^S^V' - j{^7''Yri^){i^7^^7urij) , (2-1) 

where we omit the second term of Ct and the coupling constant Gt is redefined as 
Gt = G/A. In the Dirac representation for the Dirac gamma matrices, we get 

= ° ). (2.2) 

Here, is the third component of the Pauli matrix. If the vacuum expectation 
value of V'7^7^rV non-zero value, the spin of quarks is aligned along the third 
axis which leads to the magnetization due to the quark magnetic moment. Thus, 
it is enough to investigate a possibility of spin alignment along the third axis. The 
generating functional of the Green function is expressed in the form of the Feynman 
path integral as 



(2-3) 



Following the standard technique of the auxiliary field method, we introduce the 
auxiliary field F^, which corresponds to the ferromagnetic condensate. First, the 
functional Gaussian integration. 



1 



J VFk exp -'-j d^x {Fk + G{i,E^Tk^)) G"! {Fj, + G{i^E^Tki;)) 



, (2-4) 
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is inserted in the expression of the generating functional in Eq. ()2-3p . Apart from the 
overall normalization factor, we obtain 

1 



Z oc y Vi/jVipDrk exp 
= J VFk exp 

= J VFk exp 



i J d^x (^V^n'^S^V - - Fkii^E^Tk^) 



d^x 



2G 



:Fi 



1. 



+-TY [In (7% - FkUsTk) + In {-J^p^ - FkUsTk)] 



J VFk exp i j 



d^x 



2G 



:Fi 



-Jl^ln [V + {FkTkf - FkPf.j^'UsTk + FklJsp^Y 



Tk 



■ (2-5) 



Here, from the first line to the second line, we carried out the fermionic path integral 
and from the second line to the third line, the term in which the momentum 
reverses to — is added, which gives the same contribution of the term with p^. 
Then, we divide them by factor 2. Further, from the third line to the fourth line, 
we used the relation j'^p^^'^Pu = P^- From the first line, the equation for Fk gives 
Fk = —G{tljU3Tkip)- Thus, the vacuum expectation value of Fk, that is the mean 
filed, represents the quark spin alignment or quark ferromagnetic condensate. 

The last two terms in the logarithmic function on the last line in Eq. ()2-5p except 
for FkTk is recast into 



-2i{pij 



' P2l') 



(2-6) 



where we used the relation [fij, (Tj] = i2ak for the cyclic (z, j. A;) = (1, 2, 3). Thus, the 
generating functional Z is rewritten as 

?2 



Zcx J VFkexp i J d^x (^-^ + ^trln + {FkTkf - 2i{pij^ - p2l^)FkTk]^ 



(2-7) 



Next, let us calculate the logarithmic term. We can calculate it as follows: 
tr In [V + {FkTkf - 2i(pi72 - p2^^)FkTk] 

^i{Pll^ -P2l^)FkTk] . (2- 



-trln(V + (i^ferfe)2)2+trln(l 



2i 



-p2 + (FfcTfc) 



Here, by using the expansion of the logarithmic function ln(l — x) 
^x^ — ■ ■ ■ , the second term in Eq. (|2-8p is recast into 

2i 



1^2 



tr In 1 



-p' + {FkTkf 



(pi7^ -p2l^)FkTk 
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tr 



1 / {2iFkTkf 
1 / {2iFkTkf ' ' 



-{pI+pD) 



4 V(-p2 + (F,.rfe)2)2 



{-{Pi+PD? 



— tr In 
2 



_ / 2FkTk^/pi+pi 

\ -p2 + {FkTk? 



where we used the relation 



(Pi7^ -P2l'^f = -{pI+pI) , 
tr(pi7^ -m^) = 



(2-10) 



with (7^)^ = (7^)^ = —1- Thus, finahy, we get the generating functional Z from 
Ea. (f2TH with 1^ and ([23]) as 

Z « y PFfc exp i J d^x(^-^ + ^trln (^V + (i^fcTfc)^ - 2FkTk^pl+pl^ 



+ ^tr In ( V + (FkTkf + 2FkTk\/pl + pi 



(2.11) 



It should be here noted that the argument of the logarithmic function in Eq . (|2Tip 
can be expressed as 



/ + {FkTkf ± 2FkTkJpj+pl = -pI + e(±)2 ^ 



= ( + {FkTkY ± 2FkTk^pi+pi 



Here, e^^^ is identical with the energy eigenvalue of the Dirac equation for if) derived 
from the first line of Eq. ()2-5p shown in RefJlSp . 

Since the expectation values of (ri) and {T2) for u and d quark field are zero, 
we retain only F3 which we denote as F simply. Thus, hereafter, we adopt A; = 3 
and we denote FkTk as Fer where e-r is the eigenvalue of for u and d quark field 
which gives e,- = ±1 as was shown in Ref llSQ . In general, the effective action F[F] 
is defined and introduced as 

Z = exp(ir[F]) . (2-13) 
Further, the effective potential V[F] is defined and introduced here as 

F[F] 



{FkTk) ± J Pi + Pi ] +Pi. 



(2-12) 



V[F] 



Jdh 



(2-14) 
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In our case, the effective potential V[F] is already obtained in Eq. ()2-lip in the above 
calculation: 



i(27r) 



In 



-Pt+P + {FerY - 2FerJpi + pi 



+ In -p'o +p' + {FerY + '^FerJpi + Pi 



(2-15) 



Here, trace is taken with respect to color and Dirac indices, which leads to Nc x 4. 
As for the flavor, we fix = 1 and multiply by 2 because the same contribution is 
obtained for Cr = 1 and —1. Further, the part of integration is rewritten as follows: 

'^6V[F] 



VF] 



6F 



'-dF 



2G 



ANr / dF 



d'^p 

i{2TTY 



F-^/pi+pi ^F + y/p'i + pi 



Pq -t- tp 



Po ^ 



(2-16) 



In order to extend the above treatment to a finite density case, the quark chem- 
ical potential should be introduced in the starting Lagrangian density: 

C + ^iip'^ip = C + i)^i-i^tp . (2-17) 



Therefore, we replace po as po + Further, to extend to the finite temperature 
system, the Matsubara formalism may be used.!^ We replace po and po-integration 
into the following: 



d'^p 
i{2TTy 



■f{po,p) 



d^p 
(2^ 



(2-18) 



where LOn = (2?^ -|- l)TrT is the Matsubara frequency and n is an integer. Here, T 
represents temperature and /U is a quark chemical potential. 

In the finite temperature and chemical potential system, we use the replacement 
()2-18p for the above effective potential in Eq. ()2-16p . Here, the Matsubara sum is 
taken by the standard method,'2lIi which results 



d^p F ± vrf+pf _ 1 f d^p F± ^Jp\ p\ 

2 



(27r)2 



2e 



(±)2 



n 



(±) ^(±) 



(2-19) 



As a result, the effective potential at finite temperature and chemical potential is 
derived as 



dF 



d^p 
(2^ 



F-^pi + pi 



p2^JP2_ 2FyJp{ + p'i 



n\_ ^ — ^ 



F + Vp'i+p'i 



p2 + F^ + 2F^pi+pi 



(+) (+) 
n, — n_ 



(2-20) 
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where n|^^ is defined by 



exp 



(2-21) 



+ 1 



Of course, n||I^^ is identical to the Fermi distribution function. Here, e^^^ has been 
aheady defined in Eq. (|2-12p and it also appears in the denominator of the above 
integrand. 

Finally, let us give the "gap" equation which is derived by 6V[F]/6F = 0, that 

is, 



5V[F] F 



SF 



G 



+ 2N, 



(2^ 



p2 J^p2 _ 2FyJp{+p^ 

+■ 



F + ^/pi+pi 



p^ + F^ + 2F^p{ +pi 



(+) (+) 
n\ — n 



(2-22) 



From the above gap equation, the ferromagnetic condensate F is evaluated. 



§3. Symmetric quark matter at finite density and zero temperature 

In this paper, we consider symmetric quark matter at finite density and zero 

(±) 



temperature. It should be noted that n_ represents the negative energy contribu- 
tion. Here, we should replace n'^^ into — 1 if the negative energy contribution 
is omitted. It is possible to introduce the three-momentum cutoff A to take into 
account of the negative energy contribution. The result for the phase transition 
point and so on is only obtained by replacing 1/G into a renormalized coupling 
1/Gr = 1/G — A? /tt"^. Therefore, the qualitative behavior is not changed in the later 
discussion, while the vacuum polarization is important for the spin polarization due 
to other interactions especially axial vector interaction. This procedure is de- 
scribed in Appendix B. Under this replacement, the relation n*^^ — 1 = is obtained 



at zero temperature. Thus, at zero temperature, we can express n 



(±) 



0{k^ 



where 9{x) is the Heaviside step function or the theta function. In this case, the 
effective potential ()2-20p is written as 



V\F\ 



2G 



+ 2Nc / dF 



(2^ 



F-y^p'i+p'i 



p2^p2 _ 2Fy/pl+pl 



.(-) 



+ 



F+y^p'i + pi 



p2 + ^2 + 2F^p(+pI 



(3-1) 
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Without loss of generality, F is assumed to be positive since the effective potential 
is an even function with respect to F. The role of the theta function in Eq. (j3-ip is 
to give the Fermi surface, in which the condition for the allowed momentum region 
is presented. For example, ps should satisfy the relation 

-fi<P3<f^- (3-2) 



Also, the magnitude of the remaining momentum components, \/pi + P2 ' should be 
satisfied: 



^<F- Jpf +pi<fi, i.e., F - n< Jpf +pi<F + n (for e 



-f^<F + ^pi+pl<fi, i.e., 0<^pl+pl<fi-F (for eW) . (3-3) 

Here, it is understood that the Fermi surface is modified from spherical form. There- 
fore, the effective potential should be calculated in the case F < fi and F > n 
separately. 

3.1. F > fi 

In the case F > fi, from the lower relation in Eq. ()3-3p . \/pi + P2 is not allowed 
to take a finite value except for 0. Thus, for the Fermi surface is closed. As a 



result, the second term of the integrand in Eq. ()3-ip does not give any contribution 
to the effective potential. Thus, the effective potential is computed as 



(2,) 



3 



p2 + i?2 - 2FyJp{ + pj 



The quark number density pg is similarly calculated as 

= ^^^^ / (C^^^'^ - = S^'^' ('-'^ 

with Nf = 2. 

From the gap equation 6V/5F = 0, if there exists a potential minimum in the 
region with F ^ /x, then the solution of the gap equation, F^ain: 

is obtained as 

-f^min — ■ l-J-Dj 

Then, the quark number density which gives the minimum of the effective potential 
is expressed in terms of the quark chemical potential fi as 

= J^P' . (3-7) 
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3.2. F < fi 



In the case F < fi, the Fermi surface is open for e^^^ Thus, the second term 
of integrand in Eq. (|3-ip gives a contribution to the effective potential as well as the 
first term considered in the case F > fi. Taking care of the range of integration and 
writing \/pi+ P2 = P±, we obtain 



V[F] 



El 

2G 



El 
2G 



37r2 



dF 



2^2 



dp ±p ±{F - p^) In 



f^+Vl^^-{F-p^y 
F-p± 



t^-F 



+ ^ / dp j_p ±{F +p_i_) In 



F + p± 



(SF'^l^i + 2fi^) + F^^arctan 



F\^ f^ + Vf^^ - F^ 
' A ^ F 



2 



(3-8) 



It should be here noted that the above effective potential is normalized so as to 
be equal to zero when F = 0, namely V[F = 0] = 0, due to the lower limit of 
integration being 0. The effective potential derived here with the above-mentioned 
normalization corresponds to calculating V[F] —V[F = 0]. If it is necessary to get 
V[F] itself, we have to add V[F = 0] which implies the contribution of the free 
massless quark gas and diverges due to the vacuum fluctuation.'^ In Appendix A, 
V[F] itself is given so as to avoid a divergence. 
The quark number density is derived as 



= N,Nf 



27r2 
67r2 



d^p 
(2^ 



^(M-6(-)) + 0(M-e},- 



(+)■ 



J dp_ip± V f^'^ - (F - P±y + J dp±p^^J p? -{F + p_i 



\//i2 -F2(f2 + 2/x2) + SF/u^arctan 



(3-9) 



The gap equation 5V/6F = has always a trivial solution, F = 0. If there exists the 
potential minimum in F < /i except for F = 0, then the gap equation is rewritten as 



7? 

± TT 





3^2 



2F„ 



^min + /^^arctan-^ 



2_p2 

min 



(3-10) 



Then, the quark number density which gives the minimum of the effective potential 
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is expressed in terms of the quark chemical potential fj, as 




^ {Fl,^ + 2^2)^^2_^2.^ + SF^inM^arctan 



F ■ 

-'- mm 




(3-11) 



§4. Numerical result 



As is well known, the NJL model is a nonrenormalizable model. Thus, in the 
original NJL model, the momentum cutoff parameter A and the strength of the scalar 
and pseudoscalar interactions Gs are necessary and is taken so as to reproduce the 
quark mass and pion decay constant. Usually, three-momentum cutoff A = 0.631 
GeV and Gs = 0.214 fm^ = 5.514 GeV-^ are adopted.I^ In our calculation, it is not 
necessary to introduce the momentum cutoff explicitly because the quark chemical 
potential /x, which corresponds to the Fermi energy at zero temperature, plays a role 
of momentum cutoff. Thus, in this model, only one parameter is included, namely 
G, which is the strength of the tensor interaction. 

As for the parameter G, if the Fierz transformation^^' from the original NJL 
model is assumed, then the following relation is obtained: 



Then, in this case, \Gs/Gt\ = 2 is obtained. However, in this paper, we regard G 
as a free parameter. As a guide to fix the parameter G, it is noted that Gt = G/4 
should be compared with Gs- Thus, we take Gt = 5 GeV~^, namely G = 20 GeV~^, 
which is comparable with the original NJL model parameter Gs- 

In Figdl the effective potential in Eqs. ()3-8p and ()3-4p for F < and F > fi, 
respectively, is shown as a function of F with G = 20 GeV'^. From top, /i = 0.1 GeV, 



(4-1) 



V[GeV''] 



-0.001 



0.004 



0.003 



0.002 



0.001 








/x=0.5 [GeV] 



0.8 



1.0 
F [GeV] 



Fig. 1. The effective potential is shown as a function of F. From top, /x = 0.1 GeV, 0.3 GeV, /x, 
0.45 GeV and 0.5 GeV, respectively. 
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Bred [T] 
■,16 




2 4 6 8 10 p^/p^ 

Fig. 2. The strength of "reduced" magnetic field is shown as a function of ps/po- 



0.3 GeV, /icr, 0.45 GeV and 0.5 GeV, respectively, are shown. In /i > /i^r — 0.406 
GeV, the ferromagnetic condensate F exists with non-vanishing value. Thus, the 
quark spin is aligned and the magnetic field appears spontaneously. The critical 
baryon density, which corresponds to one third of the critical quark number density, 
is about 3.47 times the normal nuclear density. Here, the critical quark chemical 
potential /icr is obtained by 6'^V[F = 0]/6F'^ = 0. Thus, the critical quark number 
density or baryon density is calculated by using the relation between the quark 
chemical potential and quark number density given in Eq. ()3-9p or ()3-5p . The critical 
density is varied as the parameter G is changed. In Table HI the critical baryon 
densities are collected. 

As is mentioned in §2, non-vanishing F means that there is a spin alignment 
of quarks, which leads to the magnetic field due to the quark magnetic moment. 
Thus, it means that there is a possibility that spontaneous magnetization occurs in 
the quark matter at high baryon density. We roughly estimate the strength of the 
magnetic field. The condensate F has the dimension of energy. Here, we introduce 
and define a "reduced" magnetic field i?red as 



red 



F 



(4-2) 



where /x^ is the average quark magnetic moment in the nonrelativistic constituent 
quark model, while the quark mass may be zero at this high density: 



^J'u + fJ'd 



2m„ 



2m„ 



(4-3) 



Table 1. The critical baryon density pa and critical quark chemical potential pci- 



G / GeV-'' 


per / PO 


/icr / GeV 


15 


5.34 


0.468 


20 


3.47 


0.406 


25 


2.48 


0.363 
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with h. Then, we obtain 

where 3mq = m^v is the nucleon mass. From Eq. ()4-2p . the strength of the "reduced" 
magnetic field is estimated. It may be expected that the "reduced" magnetic field 
corresponds to the magnitude of ferromagnetization which occurs spontaneously in 
the high density quark matter. At high density, the order of the strength of "reduced" 
magnetic field is about 10^^ ~ 10^^ T as is shown in Figl2j 



§5. Summary 

Quark ferromagnetization has been reinvestigated following the stud}^^''^ per- 
formed by two of the present authors (J. da P. and C. P.) in the framework of the 
NJL type effective model of QCD. In this paper, the effective potential with respect 
to the condensate of quark spin alignment has been derived. Actually, it has been 
shown that, at high density, quark spin alignment occurs which leads to the sponta- 
neous magnetization of the quark matter in the framework of the effective potential 
approach with the auxiliary field method based on the standard field theoretical 
technique. It is expected that non-vanishing condensate at high density leads to 
the effective field strength of spontaneous magnetization. In this paper, assuming 
the quark magnetic moment as that of the nonrelativistic constituent quark model, 
the strength of "reduced" magnetic filed has been estimated where a rather large 
magnetization has been obtained. 

In this paper, having in mind the description of the interior of compact stars, 
the quark matter at zero temperature has been treated. Of course, it is interesting 
to investigate the system at finite temperature with the Fermi distribution function 

. In this paper, we neglect the vacuum polarization by setting n^'' into — 1 
simply to regularize the divergent integral. It may be interesting to study the effect 
of the vacuum polarization or of the Dirac sea for the spin polarization in quark 
matter at finite temperaturcl^ Further, at high density, the quark matter may reveal 
the color-superconducting feature. It is also interesting to examine the coexistence 
or competition between color-superconductivity and ferromagnetizatioiJ^''^ in the 
framework developed in this paper. These are future problems. 
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Appendix A 

Reinvestigation of the effective potential 



In this Appendix, we present one comment. In the definition of V[F] in Eq. (j3-8p 
with the lower limit being of integration with respect to F, this effective potential 
actually corresponds to T^[F] — y[F = 0]. If V[F] itself is necessary as the effective 
potential, then we can express it as 



V[F] = {V[F]-V[0]) + V[0] . 



(Al) 



Here, y[F]-y[0] is finite as is seen from Eq.dMD- However, V[0] = - f d'^p/{i{2TT)'^)- 
trln7^p^ is a part of the effective potential for the free massless quark gas and 
diverges due to the vacuum fluctuation. Therefore, in order to obtain the finite 
effective potential, it is necessary to replace the last term in Eq. ()A-ip . V[0], in an 
adequate way. One possibility is that V[0] is replaced to the pressure pf^cc of a free 
massless quark gasP2} 



V[0] 



-Pb, 



180 6 Vr; 127r2 \T 



Thus, we obtain a finite effective potential at T = as 



(A-2) 



Y\F] = _ _ i!L 
2G 37r2 



F2 



(3F^/i + 2^^) + F/u^arctan 



F 



fi + ^f,^- F2 

-—r In 

4 F 



(A-3) 



with Nf = 2. In this treatment, the effective potential for F > /i should be replaced 
by 



V[F] 



El 

2G 



—fi^F 
Git 



(for F > fi) . 



(A-4) 



V[ GeV" 
0.020 




11=0.1 [GeV] _ //jO.3 [GeV] 

!1=0.5 [GeV] 

F [GeV] 



ix=0.6 [GeV] 



Fig. 3. Tlie effective potential is sliown as a function of F. From top, /x = 0.1 GeV, 0.3 GeV, fia, 
0.5 GeV and 0.6 GeV, respectively. 
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piGeV*] F„i„> M plGeV '] 




In FiglSl the effective potential defined by ()A-3p and ()A-4p is shown as a function of 
F. The gap equation is not changed. Thus, the results obtained in §4 are retained. 

The effective potential itself may be regarded as the pressure p in thermody- 
namics as is mentioned in Eq. ()A-2p already, namely, 

p = -V[F] . (A-5) 

In FigUl the pressure is shown as a function of (a) the quark chemical potential 
and (b) the baryon density divided by the normal nuclear density pQ = 0.17 fm^^, 
respectively. As is seen from FigS] (a), for quark chemical potential < Hcr, the 
branch with Fmin = is realized. From to n = 0.565 GeV, partial magnetization 
with -Fmin < /U, in which the Fermi surface is close for the states with the single 
particle energy ep^'' and these states do not contribute the magnetization, is realized. 
For large quark chemical potential, full magnetization with -Fmin > ^ is realized. As 
is seen from FigjJ] (b), at low density, the normal phase with F = is realized. 
On the other hand, at high density, the quark ferromagnetic phase with i*" 7^ is 
realized. At the intermediate region between the low and the high density, partial 
magnetization is realized. 

Appendix B 

The effect of vacuum polarization 

In this paper, we neglect the effect of the vacuum polarization because we re- 
place n^'' into n^'' — 1, in which represents the negative energy contribution, 
as was mentioned in §3. It may be important to investigate its influence for the 
ferromagnetization due to the axial vector interaction because the finite quark mass 
is essential in order to realize the ferromagnetization.^' The existence of the finite 
quark mass is due to the chiral symmetry breaking which is generated by the quark- 
antiquark condensate. Thus, the effect of the vacuum polarization, which plays a 
crucial role for generating the chiral condensate, may be important. Namely, the 
Dirac sea plays an important role to realize the chiral condensate and the dynamical 
quark mass. However, the ferromagnetization occurs without quark mass in the case 
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of tensor interaction under investigation in this paper. Therefore, we simply subtract 
the negative energy contribution. It is weh known that the role of the vacuum in 
models of the NJL type is very important and should not be ignored. The contribu- 
tion E'vac = Vvac[-^, ^] of the negative energy states, regularized by a cutoff A, can 
be easily obtained in our model. If we approximate this contribution by 

^^l , 6'K,,[F = 0,A] 

^vac ^F ^ , (B-1) 

which seems reasonable, then the effect of the vacuum is the renormalization of the 
coupling constant G. 

In this Appendix, we focus on the transition point from normal to ferromagnetic 
phase. At zero temperature, = 1. Thus, we can calculate (j2-20p with = 1. 
As for the transition point, it is only necessary to calculate 6'^V[F]/6F'^\f=o. We 
easily calculate it as 

6^V[F] 



6F^ 



F=0 ^ ^ 



with Nc = 3. Thus, the result for the phase transition point is only obtained by 
replacing 1/G into a renormalized coupling l/G-r = 1/G — vl^/vr^, namely. 



6^V[F] 



Gr 



1 1 A^ 

g; = g~^- ^^'^^ 

Therefore, the qualitative behavior is not changed in the later discussion in the 
text, while the vacuum polarization is important for the spin polarization due to 
other interactions,-^ especially axial vector interaction. If we adopt the standard 
value for the three momentum cutoff A, namely A = 0.631 GeV,'^ the critical 
baryon density /Ocr and critical quark chemical potential obtained numerically, 

which are summarized in Table |TI] instead of Table |T] given in §4. It is seen that 
the contribution of the vacuum polarization to the critical point prefers the phase 
transition at rather low density. Here, the value G = 11.028 corresponds to the 
strength obtained by the Fierz transformation as is discussed in §4. Further, around 
F = 0, the effective potential itself is expressed as 

Then, the vacuum polarization effect can also be renormalized to the coupling in 
this case. 



Table II. The critical baryon density pcr and critical quark chemical potential /icr by taking account 
of the vacuum polarization. 



G /GeV 


Pcr/pO 


Mc. /GeV 


10 


4.52 


0.443 


11.028(= 2Gs) 


3.50 


0.407 


15 


1.32 


0.294 
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